Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 6/22/04 



O 
O 
(N 

;h 

< 
in 

(N 



> 

o 
o 

^: 

6 



^ 



DETERMINING DISTANCES TO CLUSTERS OF GALAXIES USING RESONANT X-RAY EMISSION LINES 

Sandor M. Molnar\ Mark Birkinshaw^ and Richard F. Mushotzky^ 

Draft version February 5, 2008 

ABSTRACT 

Bright clusters of galaxies can be seen out to cosniological distances, and thus they can be used to 
derive cosmological parameters. Although the continuum X-ray emission from the intra-cluster gas is 
optically thin, the optical depth of resonant lines of ions of heavy elements can be larger than unity. 
In this Letter we study the feasibility of deriving distances to clusters of galaxies by determining the 
spatial distribution of the intra-cluster gas from X-ray imaging and the optical depth from resonant 
emission lines (the XREL method). We solve the radiative transfer problem for line scattering in the 
hot intra-cluster gas using Monte Carlo simulations. We discuss the spatial and spectral resolutions 
needed to use the XREL method for accurate determination of distances, and hence cosmological 
parameters, and show that accurate distances will be obtained by applying this technique with the 
next generation of high resolution X-ray spectrometers. 
Subject headings: cosmological parameters - cosmology: theory - X-rays: galaxies: clusters 



1. INTRODUCTION 

Direct physical distance measurements are of 
paramount importance in cosmology, since they 
provide us with a reliable distance scale, and make it 
possible to confirm or rule out theoretical cosmological 
models. Direct distance measurements are free from 
the systematic e rrors inheren t in the cosmic distance 
ladder method IjWebbI I1999D . Clusters of galaxies, 
the known largest gravitationally bound structures in 
the Universe, have been used extensively to determine 
cosmological parameters taking advantage of the fact 
that their formation and evolution are sensitive to 
the underlying cosmology (Allen et al. 2004; Patrick 
2004; Ettori, Tozzi and Rosati 2003; Vikhlinin et al. 
2003; Viana et al. 2003; Holder, Haiman and Mohr 
2001; Haiman, Mohr and Holder, 2001; for a review see 
Rosati, Borgani and Norman 2002). Clusters of galaxies 
can also be used as "standard candles" to determine 
cosmological parameters via distance measurements in 
the same way as the supernova method (Bonamente et 
al. 2005; Carlstrom, Holder and Reese, 2002; Mason, 
Myers, and Readhead 2001; Pointecouteau et al. 2001; 
Komatsu et al. 1999; for future prospects see Molnar, 
Birkinshaw and Mushotzky, 2002; Molnar et al. 2004). 
The next generation of such programmes aim to deter- 
mine cosmological parameters to a statistical accuracy 
of a few percent, and hence are likely to be limited 
by s ystematic errors f Molnar et al.' 2004; 'Perl mutteJ 
l2003t Haiman et all |2001; Holder ct al. 2001) New 
distance measurement techniques arc valuable because 
they permit cross-checks of measurements of cosmo- 
logical parameters made using completely independent 
methods. 

In this Letter we propose to exploit a new technique, 
the X-ray emission line (XREL) method, to determine 
direct distances to any object which, potentially, has op- 
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tical depths in x-ray resonance lines on the order of 1-5. 
This method is similar to the Sunyaev-Zel'dovich - X-ray 
(SZX) method (Sunyaev and Zel'dovich, 1980; for recent 
reviews see Carlstrom, Holder and Reese, 2002; Birkin- 
shaw, 1999), and to the X-ray resonant line absorption 
method, which is based on an absorption feature in reso- 
nance lines in s pectra of bright AGNs behind the cluster 
l)Sarazinlll98"9t iKrolik and Ravmondlll988() . We discuss 
the applicability of this technique to massive clusters of 
galaxies via Monte Carlo simulation of radiative transfer 
in the Fe He-like resonant line. 

2. X-RAY EMISSION LINES 

X-ray spectral lines contain a plethora of information 
about the p hysical and chemic al conditions of the intra- 
cluster gas llMu shotzkvl I2003D. As has been p ointed 
out by iGil'fanov. Sunvaev and Churazovl l)1987l) . reso- 
nant lines of ions of heavy elements can have opti- 
cal depths larger than unity in clusters. Photons in 
these lines suffer multiple scatterin g, and the line shapes 
and int ensities change as a result IjGil'fanov et alJ[l987t 
iShigeva ma 1998), so that the line amplitudes are reduced 
at the core of the cluster and enhanced in the outer re- 
gions. Ignoring this effect can cause an un derestimate 
in the abundances based o n resonance lines IjShigevamal 
ll998HAkimoto et al.ll2Q( ^( J|). Scatterings can also polar- 



ize th e emission lines l)Sazonov. Churazov and SunvaevI 
I2OOI . 

The line shape is also affected by gas motion in clusters. 
This can re duce the optical depth of all lines by Doppler 
broadening IjChurazov et alJ2004D . In the center of Virgo 
cluster (M87), Mathews, Buote and Brighenti (2001) es- 
timated that turbulent motion with a Mach number of 
0.3 can reduce the optical depth substantially (from 2.5 
to 1 in the case of the Si line). Sakelliou et al. (2002), 
using XMM-iVewioTi observations, found no evidence for 
resonant scattering in M87. On the other hand Xu et 
al. (2001) found evidence for resonant line scattering at 
NGC 4636 in the outskirts of the Virgo cluster. In the 
Perseus cluster resonant scattering was ruled out based 
on an analysis of the Fe XXVI K lines (Gastaldello and 
Molendi 2004; Churazov et al. 2004). However, Sanders, 
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Fabian and Dunn (2005) argue that the existence of Ha 
filaments in the core of Perseus (Conselice, Gallagher and 
Wyse 2001) suggests no turbulence. Based on resonance 
Hues of Fe XXIV and Fe XXIII, Sanders et al. (2005) 
conclude that, in Perseus, resonance scattering may be 
the cause of the apparent central drop in metal abun- 
dances. Clearly, more observations are needed to clarify 
the role of turbulence in clusters. The effects of turbu- 
lence and resonant scattering can be separated, in prin- 
ciple, since scattering broadens only resonant lines with 
non-negligible optical depth, while turbulence broadens 
all lines. Mild turbulence (up to Mach numbers of a few 
tenths) may leave some lines with optical depth larger 
than unity, so that it remains possible to determine the 
optical depth from the line shapes. In such cases we could 
use the shapes of optically-thin lines as a template and 
calculate the broadenin g due to resonance line scatter ing 
(see Figures 1 and 2 in llnogamov and Sunvaevll2003j) . 

3. THE X-RAY - RESONANCE EMISION LINE METHOD 

The XREL method compares the central X-ray sur- 
face brightness of a cluster generated by thermal 
bremsstrahlung, Sxo oc n^Q Ae(Te, Zq},) L (1 -I- z)~'^, 
where rieo is the central electron number density, Ae is 
the X-ray emissivity as a function of electron tempera- 
ture and heavy element abundance, and L is a line-of- 
sight size of the cluster, with a quantity linear in rieo L. 
In the SZX method that quantity is the central SZ effect 
or a measure proportional to it. In the XREL method we 
use the optical depth from the center of the cluster to the 
observer in the chosen resonant emission line, tq oc UeoL. 
The ratio Tq /Sxo then provides a measure of L which 
can be compared with the angular size of the cluster to 
estimate the angular diameter distance. Da, at the clus- 
ter redshift z. 

The determination of D^ in this way requires an un- 
derstanding of the shape of the cluster, and this is de- 
termined from the the X-ray image. In this Letter we 
study the main features of resonant line scattering us- 
ing a spherical isothermal beta mod el for the cluster gas 
(jCavaliere and Fusco-FemiandllQTq) . A later paper will 
study scattering in atmospheres drawn from numerical 
simulations, to investigate the effect of more realistic 
3D cluster shapes and cluster evolution. Assuming an 
isothermal beta model. Da can be expressed as 



Da'x(1 + z) FiOc, 13) Ae{Te, Zab) T^ S 



XO' 



(1) 



where F{6c,(3) is a form factor that depends on the an- 
gular core radius of the cluster, 9c = Re/ Da, and the 
shape parameter j3. Once the angular diameter distance 
- redshift function, Da{z), is determined by applying 
Equation Q] for several clusters, it can be used to con- 
strain cosmological parameters IjMolnar et alJl200^ . 

If the sensitivity and the spatial resolution of the in- 
strument are high enough, multiple distances to a clus- 
ter can be obtained by studying several resonant lines 
in different regions in the image of the cluster. The self 
checking that this provides can test for systematic errors 
in the XREL method. Determination of optical depths 
using several resonance lines and/or line intensities could 
provide additional tests. 

In the case of non-spherical clusters, the XREL method 
is subject to the same limitations as the SZX method, in 



particular to biases involved in the comparison of line-of- 
sight and angular extents of clusters. Thus an alternative 
application of the XREL method might be to help in de- 
projecting cluster density and temperature distributions 
if a cluster is at a well-known distance and suspected to 
be non-spherical. 

4. DETERMINATION OF OPTICAL DEPTH 

The optical depth in a resonant emission line can, in 
principle, be determined from the line profile, from the 
spatial distribution of the line's intensity, or from the po- 
larization of the lines, and used to determine distances to 
clusters of galaxies (Sazonov et al. 2002; Sarazin 1989). 
In this Letter we focus on the determination of the optical 
depth based on the line profile. 

The most prominent resonant line in hot clusters of 
galaxies is that of He-like iron (Fe XXV) at 6.7 keV, 
which can be detected at cosmological distances. The 
optical depth of this line can be larger than unity, al- 
though there are significant variations in the estimates 
in the literature (compare, for example Krolik and Ray- 
mond 1988 and Sarazin 1989). The spatial distribution 
of the emitted photons in the rest frame of the ion is a 
mixture of an isotropic and a dipole pattern with weights 
depending on the change in the total angular momentum 
between the two energy states participating in the tran- 
sition. 

We have solved the radiative transfer problem of line 
scattering in the intra-cluster gas via Monte Carlo simu- 
lations. This problem is well suited to Monte Carlo sim- 
ulation since the optical depth in the resonant lines is not 
excessive. We determined the line profile of the 6.7 keV 
Fe line emerging from clusters of galaxies for spherically- 
symmetric isothermal beta models for a range of optical 
depths. To, from 0.1 to 5. We carried out simulations 
using different values of /3 within the range 0.6 — 0.9 sug- 
gested by cluster data. We found that the line profile 
emerging from the core of a cluster is sensitive to tq, not 
to /3 (differs < 1%, for fixed Tq), thus, for our purposes in 
this Letter, we present results for only /3 — 2/3. We also 
included the natural hue width (FWHM == 0.335 eV), 
which is about 10% of the thermal line broadening for 
hot (5 - 10 keV) clusters of galaxies. We assumed two 
temperatures 5 and 10 kev for our simulations, and we 
worked with frequency shifts normalized by the Doppler 
thermal width, Aa; = (y — vq)/ /S.vu, and lengths nor- 
malized to the cluster core radius. Since the 6.7 keV Fe 
line originates from a pure dipole transition, there is a 
correlation between photon propagation and the polar- 
ization, and so the polarization character of the emission 
lines was included in our code. Further details of the 
code can be found in lMolnar and Birkinshaw! Ijl999() . In 
the high temperature, low density, intra-cluster gas, col- 
lision times are much longer than the spontaneous emis- 
sion times. Thus, to an excellent level of approximation, 
we may assume that photons once emitted, are not de- 
stroyed, but scatter and escape from the cluster. 

In Figure n we show the emergent line profiles for dif- 
ferent optical depths, tq = 0.5, 1, 1.5, 2, 3, 4, 5 (5 is the 
broadest) for the integrated X-ray emission from the cen- 
ter of a cluster with impact parameters, d, less than half 
the core radius {d < 0.5 Re). The dashed line represents 
the thermally-broadened Gaussian line profile character- 
istic of optically-thin emission. This figure demonstrates 
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Fig. 1. — Emerging line profiles from spherical isothermal beta 
models at a range of optical depths. See text for details. 



the diffusion of photons in frequency space. As a result 
of scattering, the line broadens. Photons can escape eas- 
ily from the wings of the line where the optical depth is 
low, but the line center will be depleted. The number of 
photons at the wings of the line is further increased by 
photons scattered into the line of sight from outside the 
nominal field of view. As a result a double-peaked line 
profile arises in the cores of clusters with optical depths 
>2, which is a unique signature of multiple scatterings. 
Comparing double-peaked lines and optically thin lines 
makes it possible to separate resonant scattering and tur- 
bulence. 

5. DETERMINATION OF COSMOLOGICAL PARAMETERS 

In principle the angular diameter distance redshift 
function, Da{z), is most sensitive to Oa at a redshift of 
around unity. Due to observational constraints, the up- 
per range of the ideal redshift window will be somewhat 
less than one. Therefore, measurements of Da{z) for a 
sample of clusters of galaxies distributed in the redshift 
range of zero to about one can be expected to provide 
useful constraints on the Hubble parameter {h = iifo/100 
km s^^ Mpc^^) and the energy densities, (O^jfiA), or 
the equation of state parameter, w = p/ p. The accuracy 
to which we can determine these cosmological parameters 
depends on the accuracy of, and the degeneracies in, the 
Da{z) function (a detailed study of errors, their degen- 
eracies and their effect on the errors in the cosmological 
parameters is given in Molnar et al. 2002). 

Errors in the angular diameter distance from the 
XREL method will be dominated by uncertainties in 
F{9c,(3), Sxo, and tq, errors from the other terms will 
be negligible by comparison. 

The accuracy of tq crucially depends on the spatial and 
spectral resolution of the instrument used. The ampli- 
tude of the dip at line center of a double-peaked line is 
sensitive to tq, and so for the best accuracy this central 
dip should be spectrally resolved. An instrument with an 
angular resolution ~ 9c/2 looking at the cluster core will 
observe a line profile close to that seen by an instrument 
with high spatial resolution. Therefore, high spatial res- 
olution is not necessary to determine tq. On the other 
hand high spectral resolution is crucial. 



Fig. 2. — The number of photons necessary to obtain 5% accu- 
racy (1 (t) in TO for spherical isothermal beta models with Tx = 5 or 
10 keV using spectrometers with 2 eV (solid), 4 eV (dashed), or 
6 eV (dot-dashed) resolution. See text for details. 



We carried out Monte Carlo simulations to study the 
feasibility of accurate optical depth determination based 
on the Fe 6.7 keV resonant line. We derived the number 
of photons in the line necessary to determine tq, with 
a fixed accuracy assuming spectroscopically-determined 
intra-cluster gas temperatures Tx (assumed to be an ac- 
curate representation of the electron temperature, Te) 
of 5 and 10 keV, which is a range for massive, ob- 
servable high redshift clusters. We collected photons 
with an impact parameter less than half the core radius 
{d < 0.5 Re). We assumed an Fe abundance of 0.3 in 
units of solar abundance, and a cluster X-ray luminosity 
- temperature scaling ( Willi s et al. 2005) . Figure|3shows 
the number of photons in the line required to determine 
To to 5% accuracy (CL 68%). Solid, dashed and dash- 
dotted lines represent X-ray observations carried out us- 
ing instruments with a spectral resolution of 2, 4 and 6 
eV (FWHM). The upper (lower) lines represent results 
for Tx = 5 keV {Tx — 10 keV). The average frequency 
change per scattering is larger for high temperature clus- 
ters, thus we need fewer photons to determine accurate 
optical depths for these clusters. Clearly, high-accuracy 
optical depth determination is feasible only for to>2, 
in which case we would need a few thousand photons in 
the Fe 6.7 keV resonant line, assuming a spectral resolu- 
tion of 2 - 4 eV. A lower-resolution spectrometer of 6 eV 
FWHM would require about ten times more line photons 
than a spectrometer with an FWHM of 2 eV. 

Practical application of this method is likely to en- 
counter a 5% statistical error iri F{6c,(3), and a similar 
error in Sxo l)Molnar et al.ll200^ . In addition, assuming 
a 5% statistical error in the optical depth, if the clus- 
ter shape is well determined, the angular diameter dis- 
tance can be determined with 10% accuracy. Measuring 
the angular diameter distance to only 100 clusters in the 
redshift interval < z < 1 with a 10% statistical error 
in Da{z), and assuming an f2„i -I- f2A = 1 (spatially flat) 
prior, would allow us to determine ftm and JIa, with 16% 
and 10% (3 a) errors. Adding an SZ survev would lead to 
a mea surement of w with a 20% (3 a) error (jMolnar et alJ 
|2004() . This is comparable to the achievable accuracy 
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based on the supernova method: using 63 supernovae, 
assuming a spatially flat model and including WMAP 
constraints, Vtm, ^k and w can be determined to ab out 
30%, 30% and 15% [la) accuracy (jKnop et alJl2003(l . 

With the next generation of X-ray instruments, it 
should be possible to obtain spectra with a few thou- 
sand photons in the Fe 6.7 keV line. Even today's in- 
struments (XMM-TVewioTT.) can collect 10^ photons in 
50 ks in the Fe 6.7 keV line complex from th e core of 
Perseus cluster finoga mov and SunvaevI 12003") . Future 
missions (Constellation-X'*, XEUS^) will have ten times 
more effective area and about 4 eV spectral resolution. 
iBonamente et all 12003) used Chandra observations of 39 
clusters (with a total of 6000 ks of data) and SZ obser- 
vations to determine the Hubble constant with a 4.5% 
statistical error. Using XMS on Constellation-X, it will 
be possible to obtain about 2000 photons in the 6.7 Fe 
Ka resonant line within half a core radius with a 15 ks 
exposure time and determine the optical depth of clus- 
ters located at about z = 0.1 with a 5% error. Taking 
account of uncertainties in the cluster shape, this should 
lead to a distance error ~ 20%, so that observations of a 
sample of 30 clusters for a total exposure time of 450 ks 
should determine the Hubble constant with about a 3% 
error. 

High redshift clusters need much longer exposures to 
obtain a few thousand photons in the 6.7 Fe resonant line 
within half a core radius of the cluster: to get useful cov- 
erage to z ^ 0.8 with Constellation-X an exposure time 
of about 1500 ks would be needed per cluster. However, 
this is likely to be an overestimate of the exposure time 
— double the number of photons per cluster can be used 
if the spatial region is increased, and for many clusters it 



should be possible to include other resonant lines (for ex- 
ample from Si and S). Also, cold core clusters have much 
higher Fe abundances at the center than we assumed in 
our simulations. Massive relaxed clusters with Tx > 8 
keV such as A1689 (z=0.183), ZW3146 (z=0.291), MS 
0451-03 (z = 0.54), and MS 1054-03 (z=0.833) seem to 
be good candidates. A detailed feasibility study is needed 
to determine the best combined spatial/spectral strategy 
and the optimal redshift distribution of the cluster sam- 
ple, but this is beyond the scope of the present Letter. 



6. CONCLUSION 

In this Letter we have estimated the achievable accu- 
racy in cosmological parameters based on direct distance 
measurements to clusters of galaxies determined from 
X-ray imaging and resonant emission lines (the XREL 
method). Our feasibility study concludes that a moder- 
ate spatial resolution (about a half of the core radius, or 
30 — 60 arcsec at the relevant redshifts) and high spectral 
resolution (ideally about 2 eV) are necessary to resolve 
the central dip in the Fe XXV 6.7 keV line produced by 
optical depth effects, and thus to achieve a useful, 10% 
cosmological distance accuracy. The next generation X- 
ray satellites will have the necessary sensitivity and res- 
olution to apply the XREL method in this way, and so 
to provide a cross-check on other methods of measuring 
cosmological parameters accurately. 

We thank G. Lake, J. Stadel and Z. Haiman for com- 
ments on our paper, P. Saha for discussions, and the 
referee, S. Ettori, for suggestions which helped to clarify 
certain points of our paper. 
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